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We develop a scaling theory to describe dynamic fluctuations of a semiflex- 
ible polymer and find several distinct regimes. We performed simulations to 
characterize the longitudinal and transverse dynamics; using ensemble averag- 
ing for a range of different degrees of coarse-graining we avoid the problems of 
slow equilibration often encountered in simulations. We find that the longitu- 
dinal fluctuations of a semiflexible object scales as t 7 ^ 8 . These fluctuations are 
correlated over a length which varies as t 1//8 . Our results are pertinent to the 
interpretation of high frequency microrheology experiments in actin solutions. 

Experiments on actin filaments have awakened interest in the dynamics of semiflexible 
polymers. Although the collective behavior in the semidilute regime has attracted much 
attention |1[], the single filament problem displays surprisingly rich dynamic features absent 
in the case of "flexible" polymers. This is a consequence of the highly anisotropic nature 
of the semiflexible polymer. The static anisotropy is characterized by the exponent for 
the growth of transverse fluctuations, r±_ as a function of the filament length, L < k |, 
r\ ~ L? I k where k is the persistence length of the filament. The dynamics is commonly 
described by a Langevin equation for the transverse fluctuations [||], 

Here, f± denotes a transverse stochastic force per unit length applied to the filament, s is 
a curvilinear coordinate and hydrodynamic interactions are neglected on the ground that 
they only induce logarithmic corrections. After time t, the filament is equilibrated over a 
length h(t) ~ (/tt) 1 / 4 . This results in a scaling of the transverse dynamics first measured in 
dynamic light scattering where it is seen that the transverse motion of a monomer on a 
filament varies as || 

(M*) 2 > ~ M*)7« ~ * 3/ 7« 1/4 (2) 

The longitudinal fluctuations of a filament are even more subtle. We need to impose strict 
incompressibility on the filament, which is violated in the simplest linearized equations such 
as (|l|). This is a long standing and delicate problem in the description of worm like polymers 
0. Recently scaling in t 3 / 4 has been predicted for the longitudinal fluctuations of a 
filament, with however an amplitude smaller by a factor L/k than in we can also 
find this result with a simple scaling argument: Due to the local conservation of filament 
length the transverse fluctuations of a segment of length I result in parallel fluctuations of 
5ri ~ | J* ds(dr±/ds) 2 comparable to (5rf) ~ 1 4 /k 2 . Each section of length li(t) of the 
filament is independent so that we can add the fluctuations of m = L/li(t) segments giving 
the motion of the end 
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(5r {l (t) 2 } ~m5rf = (L/h)(lt/K 2 ) = Lt 3/i /K 5 / 4 (3) 

While applicable to the case of a filament in a macroscopic shear field |fl this result is 
problematic for local probes of the dynamics. In particular the divergence of the fluctuations 
with the filament length must surely break down at some point. A similar problem exists in 
the athermal dynamics of semiflexible filaments P,|10l where the longitudinal friction of the 



filament modifies the dynamics. It is the purpose of this letter to explore the importance 
of longitudinal dissipation and the incompressibility constraint in the dynamics of the fully 
thermalised chain. 

In the derivation leading to (^), longitudinal friction is absent and all the different seg- 
ments of length l\(t) contribute by their longitudinal fluctuations to the motion of the end. 
However, the shortening (or extension) of a segment also requires the longitudinal motion of 
its neighbors. As a consequence, longitudinal friction limits the number of segments which 
can contribute within a finite time. Typically this limitation should not be significant if the 
longitudinal diffusion of the whole polymer (diffusing as r 2 = t/L) is so fast that there is 
no hindrance to the cumulation of fluctuations. This introduces a criterion for the validity 
of Eq. (|): Lt 3//4 //t 5//4 < t/L or L < l 2 (t) = t 1//8 K 5//8 . For longer filaments only a section of 
length l 2 (t) can contribute to the fluctuation of the end point. This limitation thus amounts 
in substituting l 2 for L in (|3]): 

(5r,|(t) 2 ) r^f^/K 5 / 8 (4) 

A similar result can be found from an analysis of the linear response of a filament 
to a weak constant longitudinal force fn using the fluctuation dissipation theorem ||11|| . 
Equation (§) predicts that in the presence of a force the end drifts as 5r\\(t) ~ f\\Lt 3 ^ / k 5 ^. 
However we cannot set the whole filament into motion at once: the typical velocity scales as 
v\\ ~ 5r\\/t ~ /||L/K 5//4 t 1//4 , and the total drag Lv\\ can not be larger than the applied force. 
This requirement again reads L < l 2 (t). For a long filament the tension propagates a distance 
hit). The corresponding section along which tension has propagated is set into motion at 
a velocity of the order of vn ~ f\\/h and the end drifts a distance 5r\\(t) ~ f^t 7 ^ 8 / k 5 ^ 8 , in 
agreement with (|4j). 

Thus a surprising conclusion is that there are two dynamic length scales which play 
important roles in the dynamics of filaments. Our main result is thus that the dynamics obey 
a pair of scaling relations with different scaling arguments. In particular, the fluctuations of 
the end of a filament are expected to follow: 

/+l/8„5/8\ 

(5) 
(6) 

In the rest of this paper we present computer simulations of a model filament in order to 
check these scaling arguments. 

We have performed simulations of a semiflexible polymer in two dimensions imposing 




a constraint on the contour length using the technique described in [12]. The polymer 



is discretized with sequence of beads with positions r*j, i e 0...n, with fixed distance 
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FIG. 1. Cloud of endpoints generated by simulating N = 100 realizations of the dynamics 
starting from the illustrated initial condition. The moments of this cloud are used to extract the 
longitudinal and transverse dynamical. We see that transverse fluctuations are larger in amplitude 
than the longitudinal fluctuations and are thus characterized by a larger moment. 

b = \ fi — fi-i\ and normalized bond vectors di = (fj — f£_i)/fe. The angles Q{ characterizing 
the bond directions are coupled by simple angular springs: 

E = lj2^0 l -0 l+ i) 2 (7) 

z 1=1 

The beads move against an isotropic friction —bdfi/dt under the influence of the forces due 
to the angular springs, —dE/dfi and stochastic forces F? a : 

at OTi 

The tensions Tj play the role of Lagrange multipliers whose values are calculated at each 
time step from the condition that the bond lengths are equal to b. 

The shortest characteristic time of the model is approximately r(b) = b 4 /n, while the 
relaxation time of a chain of n — L/b segments varies as r n ~ n 4 r(b). The total simulation 
time needed to equilibrate a chain is proportional to n 5 . The equilibration of long chains 
becomes quickly impossible; previous studies have been limited to chains which are too short 
to study detailed features of the dynamics. 

We get around the problem of generating independent configurations by simulating long 
chains for a time far shorter than the equilibration time but then performing ensemble 
averages over many short runs. These simulations are useful because we can easily prepare 
fully equilibrated initial conformations for the energy eq. ([/]) by drawing bond angles 89 = 
&i — randomly from a Gaussian distribution P{59) ~ exp {~%86 2 ^j. The choice of the 
segmentation then determines a window of accessible time scales. The elementary time step 
of the integrater is 10 _2 r(6); times shorter than r(6) are affected by discretization errors. We 
generated data for times up to 10 3 r(6) with a computational effort proportional to only n 1 . 
We chose a sequence of segmentations bj = 2~ j k and studied chains of length k/8, k/4, k/2, k 
with the number of segments varying between n = 8 and n = 1024. The overlap between 
adjacent time windows (r n = 2 4 r n / 2 ) provides a convenient check on the coarse graining 
procedure and the scaling of the parameters. 

We note that imposing a constraint on the bond lengths provides access to longer simula- 
tion times than simulations done with stiff longitudinal springs; we do not need to simulate 
the corresponding fast but uninteresting modes which limit the integration step in any ex- 
plicit integration scheme. There are nevertheless complications: An extended discussion can 
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FIG. 2. Amplitude of fluctuation as a function of time for an incompressible filament (black 
symbols, L = k) and a compressible filament (grey symbols, L = re/4, modulus K = 10 7 k _1 ) 
determined from the moments Ai(t), ^(t) of the 2D clouds. The longitudinal mode scales as 
A2W = (^ r jj) ~ i 7//8 , the transverse modes obey X\(t) = (5r"j_) ~ i 3 / 4 . For a compressible filament, 
a Rouse-like scaling \2(t) = i^ r f\) ~ t 1 ^ 2 ls found for very short times. The different symbols 
correspond to different levels of coarse-graining. 



be found in |T!|. In particular, the act of passing from elastic springs to incompressible 
bonds changes the configuration space for the problem from Euclidean to Riemannian and 
the distribution function changes from exp(— E) to V A n exp(-E). A n is the determinant of 
the Jacobian describing the transformation from Cartesian to bond angle coordinates which 
can be calculated from a transfer matrix 

A m \ = (2 -cos 2 (^-# m )\ / A, 
Aj \1 J U,;-i 

starting from the initial vector (2, 1). To recover the distribution function exp(— E) for stiff 
springs from a simulation of a rigid rod model, it is necessary to add a pseudo potential [|T3 



| log (A n ) to E and include the corresponding forces in (|8|). In the simulations described 
in this paper a proper calculation of these forces is essential to ensure that we start our short 
runs from initial conformations which are properly equilibrated. 

The objective of our simulations is the characterization of the transverse and longitudinal 
motion of the chain ends. For this purpose we perform N simulations (typically iV = 
1000) starting from an identical pre-equilibrated conformation; each simulation uses an 
independent series of random forces. We then record the N coordinates of one end of the 
chain as a function of time which form an evolving two dimensional cloud in the (x, y) 
plane. The cloud of points, Fig. 1, can be characterized by two important invariants, its two 
moments of inertia \\(t) and Xzit) (Ai > A2). The evolution of the A«(t) characterizes the 
transverse and longitudinal movement of the end of the chain (indeed Aj ~ 5r 2 ). We prepare 
a total of M random realizations (typically M = 100) of the initial chain over which we can 
calculate average properties of the cloud performing a total of MN = 10 5 simulations, Fig 2. 
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FIG. 3. Scaling behavior of the two smallest moments A^^ of the 4D clouds. Plotted is 
2{3,4} = ^{3,4} K 5 / 8 i~ 7 / 8 as a function of the tension propagation length l 2 (t) = t 1 / 8 ^ 5 / 8 for L = k 
(O), k/2 (+), k/4 (□), k/8 (x). The plot demonstrates that the longitudinal fluctuations follow 
the scaling form of Eq. (||). 



For short times the evolution of the cloud is very anisotropic, and the transverse dynamics 
corresponds to the larger moment Ai(t) which scales according to eq. @. The smaller 
moment X 2 (t) characterizes the parallel motions of the filament and for short times varies 
in agreement with eq. (^). For long times with l 2 > L, a crossover to free diffusion of the 
whole filament occurs. 

In order to understand better the relative importance of internal modes and center of 
mass diffusion in contributing to Eq. ([5]) we examined the joint motion of the two end- 
points of a chain. In a series of M simulations one generates a distribution of points in 
four dimensions (4D): {(xi, x 2 , y%, 2/2)}- The cumulants of this distribution can be used via 
the fluctuation dissipation theorem to calculate the response of the chain ends to arbitrary 



combinations of end forces pj]| . We are interested in the evolution of the four moments 
{Ai > A2 > A3 > A4} of the 4D cloud which characterize the dynamics of the whole chain. 

Our picture of the propagation of tension fluctuations as introduced above suggests the 
following scenario: As long as l 2 < L, the movement of the ends are uncorrelated. The 



4D distribution factorises and reduces to a product of the 2D case discussed above, |]T5 
the two (degenerate) smaller moments A3 = A4 ~ A2 scale as t 7//8 //t 5//8 and the two larger 
moments Ai = A2 — Ai scale as t 3 / 4 /^ 1 / 4 . For longer times l 2 (t) > L the two ends see each 
other as tension propagates along the filament. This lifts the degeneracy between the two 
smaller moments, with A 4 now characterizing the end to end fluctuations so that according 
to (3) A 4 ~ Lt 3//4 /fi; 5//4 , and A 3 characterizing the longitudinal free diffusion of the chain, 
A3 ~ t/L. Figure 3 shows the moments A3 and A4 plotted normalized by t 7 ^ 8 /K 5 ^ 8 as 
functions of t 1//8 K 5//8 /L, so that they clearly follow the scaling form Eq. (5). 

Real filaments have a longitudinal compressibility. In order to study the effect and 
to compare different methods of modeling the dynamics we have repeated some of our 
simulations using harmonic potentials E b = Yh=i ~ ~ b) 2 /2b representing elastic 
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bonds with a large but finite modulus K. These simulations require an elementary time 
step of 10~ 2 -j^2 T ip) are performed without the pseudo potential. The internal longitudinal 
dynamics are over-damped and Rouse-like {d t r\\ ~ Kd 2 r\\), giving fluctuations which scale 
as 5r 2 ~ (t/K) 1 / 2 ; the tension is correlated over a distance l 3 (t) ~ (tK) l l 2 . Compressibility 
thus introduces a third dynamic length /3 in competition with those already presented. For 
very short times t < k 5 / 3 /^ 4 / 3 the compression modes dominate over the longitudinal modes 
discussed in the first part of the article as shown in Fig 2. For longer times a cross-over to 
the dynamic behavior of the constrained chain occurs. If we assume actin filaments to be 
uniform rods of Young's modulus E and radius a, k ~ Ea A while K ~ Ea 2 . We expect 
the crossover to occur for t ~ a 8 / 3 K 1//3 and at a length (a/nY^n. For actin we estimate 
t^ 1 ~ 1 MHz, indicating that on time scales larger than 10~ 6 s the exponent 7/8 should be 
observable for filaments longer than 0.1k. 

We have presented results in two dimensions where bending and torsion are decoupled. 
In three dimensions we expect the coupling between these modes produces an even richer 
behavior [HJ. In the case of actin, a filament of length L = k has a torsional mode with 



relaxation time ~ 1ms |T7| again relaxing with a diffusive, Rouse like mode together with 



a new competing length scale. This introduces more features which we consider in a future 
publication. 

Can longitudinal fluctuations of semiflexible filaments be studied experimentally? Fol- 
lowing the motion of the end of a single filament in optical microscopy may be feasible but 
could prove somewhat tricky. Some form of scattering experiments sensitive to the longi- 
tudinal motions would be preferable. Normal dynamic light scattering measurements are 
not sensitive to longitudinal motions, however one might have some hope of introducing 



optical inhomogeneity |18] in actin filaments by marking with a low concentration of dye. 
Alternatively one might hope to immuno-attach nanometric gold beads to an actin filament 
and follow their movement with either video or light scattering techniques. 

Several of experimental groups have studied the high frequency fluctuations [|TU|-|2"2"|| of 
actin solutions with beads between 0.5/mi and 5/im in radius. The results have been in- 
terpreted using a theory for the macroscopic modulus of a solution of semiflexible objects 
starting from eq. (|). We see here that no physical cut off is needed to explain the data, 
longitudinal friction leads to a dynamic length scale ^(t) which cuts off the divergence with 
L. The beads in these experiments are sufficiently small that they should sample the point 
response of the chain rather than being sources of uniform shear. We have seen that the 
response function of a single chain has a rich scaling behavior; any experiment will sample 
a mixture of Sr 2 ~ t 1 , t 3//4 , t 1//2 and t 7 ^ 8 due to the averaging over filament lengths and 
modes. Indeed the data when examined carefully show some curvature when plotted on a 
logarithmic scale. However, for dense solutions multi-chain hydrodynamic interactions are 
surely important, we do not have a detailed understanding of this regime. 

Eventually tension propagation may be an issue in the understanding of experiments 
analyzing motor-filament systems, as it may e.g. control the response of a filament to the 
"kicks" of a collection of motors. In particular the propagation of the tension in t 1//8 will 
lead to an effective low frequency filter due to propagation delays. 

The data in this article correspond to approximately 3 x 10 6 independent simulations each 
lasting an average of 2 seconds, using 3 months of CPU time on a UltraSparc workstation. 
We would like to thank Paul Chaikin, Fred Gittes, Fred Mackintosh and Jacques Prost for 
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discussions on this work. 
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